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Abstract We consider a two dimensional skip-free reflecting random walk on a non- 
negative integer quadrant. We are interested in the tail asymptotics of its stationary 
distribution, provided its existence is assumed. We derive exact tail asymptotics for 
the stationary probabilities on the coordinate axis. This refines the asymptotic results 
in the literature, and completely solves the tail asymptotic problem on the stationary 
marginal distributions in the coordinate and diagonal directions. For this, we use the 
so-called analytic function method in such a way that either generating functions or 
moment generating functions are suitably chosen. The results are exemplified by a 
two node network with simultaneous arrivals. 



1 Introduction 

We are concerned with a two dimensional reflecting random walk on a nonnegative 
integer quadrant, which is the set of two dimensional vectors such that i,j are 
nonnegative integers. We assume that it is skip free in all directions, that is, its incre- 
ments in each coordinate direction are at most one in absolute value. The boundary 
of the quadrant is partitioned into three faces, the origin and the two coordinate axes 
in the quadrant. We assume that the transition probabilities of this random walk is 
homogeneous on each boundary face, but they may change on different faces or the 
interior of the quadrant, that is, inside of the boundary. 

This reflecting random walk is referred to as a double quasi-birth-and-death pro- 
cess, a double QBD process for short, in lfT8ll . This process can be used to describe 
a two node queueing network under various setting such as server collaboration and 
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simultaneous arrivals and departures, and its stationary distribution is important for 
the performance evaluation of such a network model. The existence of the stationary 
distribution, that is, stability, is nicely characterized, but the stationary distribution 
is hard to analytically get except for some special cases. Because of this as well as 
its own importance, research interest has been directed to its tail asymptotics. 

Until now, the tail asymptotics for the double QBD have been obtained in terms 
of its modelling primitives under the most general setting in Miyazawa lfT8l . while 
less explicit results have been obtained for more general two dimensional reflecting 
random in Borovkov and Mogul' skii (2). Foley and McDonald ifTOlfTTI studied the 
double QBD under some limitations. Recently, Kobayashi and Miyazawa [ 13 1 mod- 
ified the double QBD process in such a way that upward jumps may be unbounded, 
and studied its tail asymptotics. This process, called a double M/G/l type, includes 
the double QBD process as a special case. For special cases such as tandem and pri- 
ority queues, the tail asymptotics have been recently investigated in Guillemin and 
Leeuwaarden lfl2l and Li and Zhao fl4l[T5l . Recently, Li and Zhao [ 16] challenged 
the general double QBD (see Additional note at the end of this section). 

The tail asymptotic problems have also been studied for a semi-martingale re- 
flecting Brownian motion, SRBM for short, which is a continuous time and state 
counterpart of a reflecting random walk. For the two dimensional SRBM, the rate 
function for large deviations has been obtained under a certain extra assumption 
in Avram Dai and Hasenbein [1|. Dai and Miyazawa [3| derived more complete 
answers but for the stationary marginal distributions. 

Thus, we now have many studies on the tail asymptotics for two dimensional 
reflecting and related processes (see, e.g., |fl9l for survey). Nevertheless, there still 
remain many problems unsolved even for the double QBD. The exact tail asymp- 
totics of the stationary marginal distributions in the coordinate directions are one of 
such problems. Here, a sequence of nonnegative numbers {p(n);n = 0,1, 2} is said 
to have exact tail asymptotic {h(n)\n = 0,1,...} if their ratio p(n)/h(n) converges 
to a positive constant as n goes to infinity. We also write this asymptotic as 

p(n) ~ h[n) 

We will find h(n) = n K aT n or n K (l + b(—l) n )a~ n with constants K = 1, 
a > 1 and \b\ < 1 for the marginal distributions (also for the stationary probabilities 
on the boundaries). 

We aim to completely solve the exact tail asymptotics of the stationary marginal 
distributions in the coordinate and diagonal directions, provided the stationary dis- 
tribution exists. It is known that the tail asymptotics of the stationary probabilities 
on each coordinate axis are a key for them (e.g., see [19]). These asymptotics have 
been studied in El 1 181 . They used Markov additive processes generated by remov- 
ing one of the boundary face which is not the origin, and related their asymptotics. 
However, there are some limitations in that approach. 

In this paper, we revisit the double QBD process using a different approach, 
recently developed in |3j [T3] |2T) . This approach is purely analytic, and called an 
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analytic function method. It is closely related to the kernel method used in lfl2l fl4l 
[T51 . Their details and related topics are reviewed in |[T9l . 

The analytic function method in (3] Q~3j |2TI only uses moment generating func- 
tions because they have nice analytic properties including convexity. However, a 
generating function is more convenient for a distribution on integers because they 
are polynomials. Thus, generating functions have been used in the kernel method. 

In this paper, we use both of generating functions and moment generating func- 
tions. We first consider the convergence domain of the moment generating function 
of the stationary distribution, which is two dimensional. This part mainly refers to 
recent results due to lfT3ll . Once the domain is obtained, we switch from moment 
generating function to generating function, and consider analytic behaviors around 
its dominant singular points. A key is the so called kernel function. We derive in- 
equalities for it (see Lemma[8]l, adapting the idea used in |3|. This is a crucial step 
in the present approach, which enables us to apply analytic extensions not using the 
Riemann surface which has been typically used in the kernel method. We then apply 
the inversion technique for generating functions, and derive the exact tail asymp- 
totics of the stationary tail probabilities on the coordinate axes. 

The asymptotic results are exemplified by a two node queueing network with 
simultaneous arrivals. This model is an extension of a two parallel queues with 
simultaneous arrivals. For the latter, the tail asymptotics of its stationary distribution 
in the coordinate directions are obtained in |8, 9|. We modify this model in such a 
way that a customer who has completed service may be routed to another queue with 
a given probability. Thus, our model is more like a Jackson network, but it does not 
have a product form stationary distribution because of simultaneous arrivals. We 
will discuss how we can see the tail asymptotics from the modeling primitives. 

This paper is made up by seven sections. In Section [2] we introduce the double 
QBD process, and summarize existence results using moment generating functions. 
Section [3] considers the generating functions for the stationary probabilities on the 
coordinate axes. Analytic behaviors around their dominant singular points are stud- 
ied. We then apply the inversion technique and derive exact asymptotics in Sections 
[4] and [5] The example for simultaneous arrivals is considered in Section [6] We dis- 
cuss some remaining problems in Section[7] 

(Additional note) After the first submission of this paper, we have known that Li 
and Zhao [16| studied the same exact tail asymptotic problem, including the case 
that the tail asymptotics is periodic. This periodic case was lacked in our original 
submission, and was added in the present paper. Thus, we benefited by them. How- 
ever, our approach is different from theirs although both uses analytic functions and 
its asymptotic inversions. Namely, the crucial step in [16] is analytic extensions on 
a Riemann surface studied in [5|, while we use the convergence domain obtained 
in lfl3l and the key lemma. Another difference is sorting tail asymptotic results. 
Their presentation is purely analytic while we use the geometrical classifications of 
C2 [13 (see also EI). 
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2 Double QBD process and convergence domain 

The double QBD process was introduced and studied in ifTHIl . We here briefly intro- 
duce it, and present results on the tail asymptotics of its stationary distribution. We 
will use following set of numbers. 

Z = the set of all integers, Z + = {j £ Z;j > 0}, 

U = {(ij)eZ 2 UJ = 0,1,-1}, 

K = the set of all real numbers, M + = {x G M.;x > 0}, 
C = the set of all complex numbers. 

Let S = which is a state space for the double QBD process. Define the boundary 
faces of S as 

<9S = {(0,0)}, dSi={(i,0)eZ 2 + ;i>l}, dS 2 = {(0,0 € %\\ i > 1}. 

Let dS = L)j =0 dSi and S + = S \ dS. We refer to dS and S + as the boundary and 
interior of S, respectively. 

Let {Yf : ;£ = 0, 1, . . .} be a skip free random walk on 1? . That is, its increments 
X^ + ' = Y( — Y^_i take values in U, and are independent and identically distributed. 
By XW, we simply denote a random vector which has the same distribution as 
X^ + '. Define a discrete time Markov chain {L^} with state space S by the transition 
probabilities: 

where X^' is a random vector taking values in { (i'i , 1*2) € U;/3_^ > 0} for k = 1,2 
and in , 12) € U; i\ , 12 > 0} for k = 0. Hence, we can write as 

U+x=U+ I xfi(L^e^), £ = 0,1,2,.... (l) 

(1=0,1,2,+ 

where l(-) is the indicator function of the statement and xf' has the same 
distribution as that of X^' for each k = 0,1,2, +, and is independent everything 
else. 

Thus, {L^} is a skip free reflecting random walk on the nonnegative integer quad- 
rant S, which is called a double QBD process because its QBD transition structure 
is unchanged when level and background state are exchanged. 

We denote the moment generating functions of X^ by 7^, that is, for 9 = 

(9 u 9 2 )eR 2 , 

7k {9)=E{e^- x(k \ k = 0,1,2,+, 
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where (a,b) = aibi +02^2 for a = (a\ ,a%) and b = (61,^2). As usual, R 2 is consid- 
ered to be a metric space with Euclidean norm ||a|| = yj (a, a). In particular, a vector 
c is called a directional vector if ||c|| = 1. In this paper, we assume that 

(i) The random walk {Ye } is irreducible. 

(ii) The reflecting process {Lf } is irreducible and aperiodic. 

(iii) Either + or£(x[ +) ) ^ for X« = {x[ +) ,X^ +) ). 

Remark 1. If E(x[ + ' i ) = Eft^) — 0, then it is known that the stationary distribu- 
tion of {L^} cannot have a light tail, that is, cannot geometrically (or exponentially) 
decay in all directions . See (6 1 and Remark 3.1 of (13). Thus, (iii) is not a restrictive 
assumption for considering the light tail. 

Under these assumptions, tractable conditions are obtained for the existence of 
the stationary distribution in the book [6]. They are recently corrected in [13|. We 
refer to this corrected version below. 

Lemma 1 (Lemma 2.1 of lfl3) ). Assume condition (i)-(iii), and let 

m=(E(x[ +) ),E(4 +) ))> 
mM = (E(4\-E(xW)), 

m® = (-E(X®),E(X®)). 

Then, the reflecting random walk {L(} has the stationary distribution if and only if 
either one of the following three conditions hold (see IT3l ). 

mi < 0,m 2 < 0, (n^m^) < 0, (m,m^ 3 ) < 0, (2) 
m\ > 0,ni2 < 0, (m,m^) < 0. In addition, m^ < is needed if mf 1 = 0. (3) 
mi < Q,ni2 > 0, (m,m^') < 0, In addition, Wj < is needed if m^p — 0. (4) 

Throughout the paper, we also assume this stability condition. That is, 

(iv) Either one of |2|, ^ or Q holds. 

In addition to the conditions (i)-(iv), we will use the following conditions to 
distinguish some periodical nature of the tail asymptotics. 

(v-a) P(xWe{(l,l),(-l,l),(0 1 0),(l 1 -l),(-l,-l)})<L 
(v-b) P(X( 1 )e{(l,l),(0,0) I (-l,l)})<l. 

(v-c) p(x( 2 >e{(M) J (o,0),(i,-i)})<i. 

These conditions are said to be non-arithmetic in the interior and boundary faces 
1,2, respectively, while the conditions that they do not hold are called arithmetic. 
The remark below explains why they are so called. 
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Remark 2. To see the meaning of these conditions, let us consider random walk 
{Y(} on Z 2 . We can view this random walk as a Markov additive process in the k-th 
coordinate direction if we consider the k-th entry of as an additive component and 
the other entry as a background state (k = 1,2). Then, the condition (v-a) is exactly 
the non-arithmetic condition of this Markov additive process in each coordinate 
direction (see [20 1 for the definition of the period of a Markov additive process). 
It is notable that, for the random walk {Y^}, if the Markov additive process in one 
direction is non-arithmetic, then the one in the other direction is also non-arithmetic. 

We can give similar interpretations for (v-b) and (v-c). Namely, for each k= 1,2, 
consider the random walk with increments subject to the same distribution as X". 
This random walk is also viewed as a Markov additive process with an additive com- 
ponent in the k-th coordinate direction. Then, (v-b) and (v-c) are the non-arithmetic 
condition of this Markov additive process for k — 1,2, respectively. 

Remark 3. These conditions were recently studied in [ 16 1. They called a probability 
distribution on U = = —1,0, 1} to be X-shaped if its support is included 

in 

{(1,1), (-1,1), (0,0), (1,-1), (-1,-1)}. 

Thus, the conditions (v-a), (v-b) and (v-c) are for XW, XW and X (2 \ respectively, 
not to be X-shaped. 

We denote the stationary distribution of {Lf,£ = 0,1,...} by v, and let L be a 
random vector subject to v. Then, it follows from ([T]i that 

L~L + £ X«1(LGS*), (5) 

fc=0,l,2,+ 

where "~" stands for the equality in distribution. We introduce four moment gener- 
ating functions concerning v. For 9 £ R 2 , 

<p + (0)=£(^ L >l(LeS + )), 
< Pk (6 k )=E(e e * L n(L£S k )), £=1,2. 

Then, from |5]) and the fact that 

2 

< p(0) = <p + (0) + £<p fc (0 t ) + v(O), 

k=l 

we can easily derive the following stationary equation. 

(i-r+(0))(p + (0) + (i-7i(0))(Pi(0i) 

+(1 - 7 2 (0))<P2(0 2 ) + (1 - *>(0))V(O) = 0, (6) 
as long as (p(9) is finite. Clearly, this finiteness holds for 9 < 0. 
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To find the maximal region for |6]) to be valid, we define the convergence domain 
of <p as 

^ = the interior of {0 e M 2 ; 9(0) < °o}. 

This domain is obtained by Kobayashi and Miyazawa fl3l . To present this result, 
we introduce notations. 

From ([6|, we can see that the curves 1 — Jk{9) = for k = +, 1,2 are keys for 
9(0) to be finite. Thus, we let 

r l = {eei 2 ;K (e)<i} ! 

«9r, = {0eM 2 ;%(0) = l}, *=1,2,+. 

We denote the closure of fj. by IV Since jk is a convex function, Ijj and Tj are con- 
vex sets. Furthermore, condition (i) implies that r + is bounded, that is, it is included 
in a ball in M 2 . Let 

9 {k ^ = arg eeR2 sup{9 k ; 9 e T+ nl*}, * = 1,2, 
e ( *' mm) =arg eGR2 inf{0,;0Gr+}, 

0^ max )=arg eGK2 sup{0,;0er + }. 

These extreme points play key roles in obtaining the convergence domain. It is no- 
table that is not the zero vector because the stability condition (iv) implies 
that, for each k = 1,2, 1+ Hi* contains = (0i, ©2) such that 0j. > (see Lemma 
2.2 of HI). 

We further need the following points. 

fl (*,n _ / o (k ' T) , %(e (Lmax) )>i, , 19 

According to Miyazawa [ 18 1 (see also [3 1), we classify the model into the following 
three categories. 

Category I 0« 2 ' r) < 9^ and 9^ ] < 9^\ 

Category H 0p' F) < 1 (1,r) and 2 (i r) > 2 (2 ' r) , 

Categoiy HI 9^ T) > 9 { l 1T) and 2 (1 ' r) < 9^ F) . 

(2 r) fi -T) fir) (2 r) 
Note that it is impossible to have 0| ' > 0j and 2 > 2 at once because 

9[ 2T) > 0! (i r) and the convexity of T + imply that 2 (l r) < 2 (2,r) (see Section 4 

of EH). We further note that 9^' F) > 2 (2 ' F) can be replaced by 2 (1,r) = 2 (2 r) in 

(2 r) "(if) (2 r) (1 r) 

Category n. Similarly, 0j > 0j ' can be replaced by 0j — 0j ' in Category 

m. 

Define the vector T as 
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{(d[ l ' r) , ftf ' r) ), for category I, 

(^(0f r) ),ef r) ), for category H, 
(9^,^(6^)), for category ffl , 

where ^ k {dj,-k) = sup{0jt; (9\ , 82) € 1+}. This definition of T shows that categories 
I, II and IE are convenient. 

We are now ready to present results on the convergence domain £? and the tail 
asymptotics obtained by Kobayashi and Miyazawa lfl31 . As we mentioned in Sec- 
tion [TJ they are obtained for the more general reflecting random walk. Thus, some 
of their conditions automatically hold for the double QBD process. 

Lemma 2 (Theorem 3.1 of lH3l). 

9 = {9 e R 2 ; 9 < x and 30' G T+ such that 9 < 9'}. (7) 




Theorem 1 (Theorem 4.2 of 1 13 1). Under conditions (i)-(iv), we have, for k — 1.2, 

lim -\ogP{L k >nM-k = 0) = -x kl (8) 

n— y°° n 

and, for any directional vector c > 0, 

lim -logP((c,L) >x) = -a c , (9) 

where we recall that a c — sup{x > 0;xc G 3>}. Furthermore, if y(ct c c) = 1 and if 
7i:(a c c) 7^ 1 and a c c k 7^ djZ° for k = 1,2, then we have the following exact asymp- 
totics. 

lim e acT P((c,L) >x)= b c . (10) 
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In this paper, we aim to refine these asymptotics to be exact when c is either 
(1,0), (0, 1) or (1, 1). Recall that a sequence of nonnegative number {p(n);n £ Z + } 
is said to have the exact asymptotic (1 +b(— l)")n~ K a~" for constants K and a > 1 
if there exist real number be [—1,1] and a positive constant c such that 

Yim(\+b{-\)")n K a n p{n)=c. (11) 

n— »oo 

We note that this asymptotic is equivalent to 

lim (1 +b{-\) n )n K a" TpU) = c' (12) 

for some d > 0. Thus, there is no difference on the exact asymptotic between > 
n) and P(Lk =n). In what follows, we are mainly concerned with the latter type of 
exact asymptotics. 



3 Analytic function method 

Our basic idea for deriving exact asymptotics is to adapt the method used in 
which extends the moment generating functions to complex variable analytic 
functions, and gets the exact tail asymptotics from analytic behavior around their 
singular points. A similar method is called a kernel method in some literature 
|[T2l[T4l[T5l[T6ll . We here call it an analytic function method because our approach 
heavily use the convergence domain which is not the case for the kernel method. 
See |[T9ll for more details. 

There is one problem in adapting the method of [ 3 ] because the moment gen- 
erating functions Jfc(0) are not polynomials, while the corresponding functions of 
SRBM are polynomial. If they are not polynomials, the analytic function approach 
is hard to apply. This problem is resolved if we use generating functions instead of 
moment generating functions. We here thanks for the skip free assumption. 



3.1 Convergence domain of a generating function 

Let us convert results on moment generating functions to those on generating func- 
tion, using a mapping from z = (zi,Z2) € C to g(z) = (e Zl ,e Z2 ) 6 C. In particular, 
for £ K 2 , g(0) € (R° ) 2 , where M° = (0,»). We use the following notations for 
£=1,2. 



10 M. Kobayashi and M. Miyazawa 

(uf min) ,4"- mm) ) = g^-™")), ( M f ,max) ,4* ,max) ) = g(e (A -- max) ), 
(«S w ,4 M) )=«(e (W ), («^,«4* ,n )=«(fl (W )- 

(Tl,T 2 )=g(T) ! 

We now transfer the results on the moment generating functions in Section[2]to those 
on the generating functions. For this, we define 

l = {g(e)eM 2 +; e€f}, 

r k = {g(O)eR 2 + ;0er k }, k=l,2,+- 
Define the following generating functions. For k = 0, 1 , 2, +, 

y(k) y(l) 

%(z)=E$ zf ), Z E( Z|)Z2 )eC 2 , 

which exists except for z\ — or Z2 = 0. Similarly, 

<p(z)=E(z L 1 l Z?), 
? + (z)=£(^^l(L6S + )), 

ft(z k )=£(z^l(LGS'i t )). fe = 1 ' 2 - 

as long as they exist. 

Obviously, these generating functions are obtained from the corresponding mo- 
ment generating functions using the inverse mapping g _1 . 

%0) = %(logzi,logz 2 ), £ = 0,1,2,+, 
<p(z) = <p(logzi,logz 2 ), 
<M Z ) = <P+(logzi,logz 2 ), 
Pk(z) = <Pft(logz), £=1,2, 

Then, the stationary equation |6]l can be written as 

(l-?+(*))0+(*) + (l-MsO)fr(zi) 

+(1 - fc(z))fc(z 2 ) + (1 - fo(z))v(0) = 0. (13) 

It is easy to see that 

A={ueR5.;u>0,ft(u)<l}, k=l,2,+, 
fE{uel 2 + ;n>fl,9(u) <°c}, 

It is notable that these sets may not be convex because two dimensional generat- 
ing functions may not be convex (see Figure [2]). Nevertheless, they still have nice 
properties because the generating functions are polynomials with nonnegative coef- 
ficients. To make this specific, we introduce the following terminology. 
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Fig. 2 The examples of S> and the corresponding 9, which may not be convex, where 
(p2l,P0l,Pli) = (0.1,0.1,0.7), (p 2 o, Poo, Pio) = (1.5,0.5,0.5), (p22,P02,Pi2) = (2,3,1) forp, ; = 
P(X<+) = (/,;)). 
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Definition 1. A subset A of Mr is said to be nonnegative-directed (or coordinate- 
directed) convex if Ax + (1 — A)y £ A for any number A G [0, 1] and any x,y£A 
such that y — x > (or y — x in either one of the coordinate axes, respectively). 

We then immediately have the following facts. 

Lemma 3. !& is nonnegative-directed convex, and is coordinate-directed convex 
for for £ = +,0,1,2. 



13 1 is valid for z G C satisfying |z| G ^ because |<p(z)| < <p(|z|). 



Note that 
Furthermore, 

{z g C 2 ; |z| el} = {ze c 2 ;^^ 1108 !-- 1 !^ 21081 - 72 !) < «.} 
= {g(log|zi|+/argz 1 ,log|z 2 |+/argZ2);zeC 2 ,(log|zi|,log|z 2 |) G $>} 
= g({zeC 2 ,(5Rz 1 ,9iz 2 )e^}), 

where |z| = (|zi|, |z2|)- Hence, the domain is well transferred to We will work 
on Q} for finding the analytic behaviors of (pi(z) and <p 2 (z) around their dominant 
singular points. This is different from the kernel method, which directly works on 
the set of complex vectors z satisfying y+(z) = 1, and applies deeper complex anal- 
ysis such as analytic extension on a Riemann surface (e.g., see Q). We avoid it 
using the domain 3. 



3.2 A key function for analytic extension 



Once the domain @ is obtained, the next step is to study analytic behaviors of the 
generating function (fa for k = 1,2. For this, we use a relation between them by 
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letting y + (z) — 1 = in the stationary equation ( 13 1, which removes <p+ (z) . For this, 
let us consider the solution u 2 > of y + (u\ , u 2 ) = 1 for each fixed u\ > 0. Since this 
equation is quadratic concerning u 2 an d C (R°)?_, it has two positive solutions 
for each u\ satisfying 

(l,min) . . „(l,max) 

Denote these solutions by C 2 ( M i) an d C 2 ( u i) sucn mat C 2 ( M i) — C 2 ( u i)- Similarly, 
Cj("2) an d Ci( M 2) are defined for m 2 satisfying 

(2,min) . . „(2,max) 

One can see these facts also applying the mapping g to the convex bounded set 5? 
(see Lemma[3]l. 

We now adapt the arguments in Q. For this, we first examine the function 
Let ~ 

p* k {u)=E{u x ^l{x( +) =*)), 

Pt*(M) =£( M X 2 +) l(X, (+) =jt)), jfc= 0,1,-1. 

Then, f+( u i> u 2) = 1 can be written as 

(hi) -«2(1 — P*o(mi))+/>*-i(mi) = 0. (14) 

Hence, we have, for m 6 [m^ mm \ui max )], 

r / \ l-p*o(u)-y/D 2 (u) 
- 2 2p*i{u) 

where 

D 2 (u) = (1 - p*o{u)) 2 -4p t i(u)p t .-i(u) > 0. 

Since £) 2 (M ( 1 1 ' min) ) = D 2 (u[ Umax] ) = and u 2 D 2 (u) is a polynomial with order 4 at 
most and order 2 at least by condition (i), u 2 D 2 (u) can be factorized as 

2 n / \ / (l.min)w M.max) \, / \ 

u D 2 [u) = (u — u\ )\ u \ — u)n 2 {u), 

where %(w) 7^ for « e (m! 1 '™',!*' 1 ), This fact can be verified by the mapping 
gfromr+ toT+. 

To get tail asymptotics, we will use analytic functions. So far, we like to analyt- 
ically extend the function £ 2 from the real interval to a sufficiently large region in 
the complex plane C. For this, we prepare a series of lemmas. We first refer to the 
following fact. 
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Lemma 4 (Lemma 2.3.8 of 0). All the solutions of z 2 L>2 (z) = for z € C are real 
numbers. 

In the light of the above arguments, this lemma immediately leads to the follow- 
ing fact. 

Lemma 5. z 2 Z?2(z) = for z <E C has no solution in the region such that \z\ € 

, (l,min) (l,max)\ 
\ u l i u \ )■ 

We will also use the following two lemmas, which show how the periodic nature 
of the random walk {Y^} is related to the branch points (see Remark [2] for the 
periodic nature). They are proved in Appendices A and B, respectively. 

Lemma 6. The equation: 

D 2 (z)=0, |z|=4 1,max) ,zGC, (16) 

has only one solution z = Mj l max ' if and only if (v-a) holds. Otherwise, it has two 
solutions z = ±«[ and u 2 D2(u) is an even function. 

Lemma 7. For each fixed x,y > 0, we have 

f + (x,y) = l, f+(-x,-y) = l, (17) 
if and only if (v-a) does not hold. 

Remark 4. Lemma[6]is essentially the same as Remark 3. 1 of ifTBI , which is obtained 
as a corollary of their Lemma 3. 1, which is immediate from Lemmas 2.3.8 of 0. 

By Lemmas ^ and j^J C 2 ( M ) on {u[ l ' mm \ i/ 1 1 ' max ') is extendable as an analytic 
function of a complex variable to the region %(«[ 1,mm ' , «[ 1 ' max ' > ), where 

%(a,b) = {ze C;z £ (-°°,a\ U [&,«»)}, a,b € E, 

and has a single branch point Mj l max ' on \z\ = M| l max ' if (v-a) hold, and two branch 

points ±i/ 1 l max ' there otherwise by Lemmas^ and [5] Both branch points have order 
two. We denote this extended analytic function by~£ 2 (z). That is, we use the same 
notation for an analytically extended function. We identify it by its argument. The 
following lemma is a key for our arguments. The idea of this lemma is similar to 
Lemma 6.3 of [3 1, but its proof is entirely different from that lemma. 

Lemma 8. (a) £ 2 of is analytically extended on %(u[ l ' mm \u[ l ' max ^). 
(b) Forz G C satisfying z\ € (i/ 1 1:mm \i<i 1 ' max '], 

\Uz)\<U\z\)<u^\ (18) 
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where the second inequality is strict if \z\ < u^ ,7mx \ 

(c) If either nffl — or (v-b) holds, then 

7i(z,C 2 (z)) = l, \z\ = u^\ (19) 

(1 r) 

has no solution other than z = u\ . 

(d) The equation |l9| ) has two solutions z — ±Wj 1,r ^ if and only if neither = 0, 
(v-a) nor (v-b) holds. 

Proof. We have already proved (a). Thus, we only need to prove (b), (c) and (d). 
We first prove (b). For this, it is sufficient to prove (18i for \z\ < i/ 1 1 ' max ' ) by the 

continuity of C 2 (k) f° r kl — Hj ' ma "^ at z = u[ l ' maK \ Substituting complex numbers 
Zi and Z2 into u\ and U2 of ( [14) , we have 

Z2P*\{Z\)+Z2P*0{Z\)+P*-\{Z\) =Z2- (20) 

Obviously, this equation has the following solutions for each fixed z\ such that |zi | G 

/ (l.min) (l.maxK 

Z 2 = C 2 (*l), C 2 (Zl)- (2D 

We next take the absolute values of both sides of ( [20) , then 

|Z2| 2 p*l(kl|) + k2b*0(|zi|)+P*-l(|zi|)> |Z2|. 

Thus, we get 

k 2 |(y+(ki|,k 2 |)-i)>o. 

By the definitions of C 7 (|zi |) and C 2 (ki I)' m i s inequality can be written as 

(k2|-C 2 (ki|))(k2|-C 2 (kil)) = k2|(y + (ki|,k 2 |)-i)>o. 

Hence, C 2 (|zi|) < C 2 (ki|) implies 

|z 2 |<C 2 (ki|) or C 2 (kil)<k2|. (22) 
By (21 1, we can substitute z 2 = C 2 (^i) into (22 1, and get 

IC 2 (zi)l<C 2 (ki|) or C 2 (kil)<IC 2 (zi)l, \zi\€(u^\u^). (23) 

Thus, (b) is obtained if we show that C 2 (ki I) ^ IC 2 (zi)| is impossible. Suppose 
the contrary of this, that is, there is a z\^ such that 
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c 2 (i4 0) i)<ic 2 (4 0) )i, 



JO), 



e (m 



(l,min) (l,max) 
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(24) 



Since | C, (z) | is continuous and converges to £ 



,(0), 



as z goes to |z[°^ | on the path 



that |z| = Izi '!, there must be a zi 1-1 such that |z 



(i) 



|zS 0) l and 



c 2 (i4 1) i)<ic 2 (4 1) )i<c 2 (k ( i 1) i) 



Since |z 



,(0), 



e (u 



(l,min) (l.max)- 



, this contradicts (23 1, which proves (b). 



We next prove (c). Let 

pU(z)=E(z x r'l(X^=k)), A=0.1. 



First, assume that mi = 0. This implies p^' (z) = 0, and therefore ( 19 1 is reduced to 

pty (z) = 1. Hence, its solution is z = 1 or z = pf!jn/ Mo ^ if p\xj~£ (otherwise, 
z = 1 is only the solution). Both are nonnegative numbers, and therefore ( 19 1 has no 
solution z such that 



(V) 



(25) 



We next assume that 7^ 0, which implies p'j'(z) ^ 0. Since (|19|) can be 



written as 



C 2 (z)lfi ) (z)+J>2(*) = l 



(26) 



and 1 < \w\ + 1 1 — w\ for any w E C, we have 

IC,(z)l 



> i-blo ) WI > i- 



= C,(W)- (27) 



If ( |25| l holds, then both sides of this inequality are identical if and only if (v-b) does 
not hold. Hence, if (v-b) holds, then |£,(z)| > C 2 (kl)> an( ^ therefore (19l has no 
solution satisfying d25]l because of (fT8|). 

(2) 

We finally prove (d). For this, we assume that both of m\ — and (v-b) do not 



(i) 



hold. In this case, p M 



.Ji) (i) 
' "(-1)0 -Pio 



0, so it follows from ( 26 1 that 



Z7 (1) WV 
P-n + Pn z 



Hence, if (|25|l holds, then we must have z 

^2 



( 1 t) 

above equation, we also have u\ ' ) 



because of (|18j and (|27j). By the 
"C<j( M i )■ Hence, we need to check 



whether (— u^ ,T \ — ^(Mj 1 '^)) be the solution of y + {x,y) — 1. By Lemma |tJ z 
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Table 1 The solutions of 




where Q, X and — indicate "yes", "no" and '"irrelevant". 



Non-arithmetic: (v-a) 





X 


X 


X 


X 


Non-arithmetic: (v-b) 










X 


X 


m ( 2 l) = 







X 





X 


The solutions of ( 


16 




(l,max) 


, (l.max) 
±U\ ' 


The solutions of ( 


19 




u\ 


(l,r) 





—u\ is the solution of |l9| if and only if (v-a) does not hold. Combining this with 
(b) and (c) completes the proof of (d). □ 

For convenience of later reference, we summarize the results in (c) and (d) of 
Lemma [8] in Table Q] Similar results can be obtained in the direction of the 2nd 
axes using (v-b) and Wj = instead of (v-a) and rti^ — 0. Since the results are 
symmetric, we omit them. We remark that Li and Zhao ifTBI have not considered the 
cases m[ 2 ' = and trip = 0, which seems to be overlooked. 



3.3 Nature of the dominant singularity 

We consider complex variable functions (p\ (zi) and ^2(^2)- Recall that 

4K«) = 4>+(*) + fr(zi) + fcfez)+v(0). (28) 

Obviously, <p(z) is analytic for z G C 2 such that (|zi|, |Z2|) € an d singular on the 
boundary of §f. This implies that <Pi(z;) is analytic for |z,-| < f; and has a point on the 



circle \z\ = fj. This is easily seen from (28 1 with Zj = for j = 3 — i. Furthermore 



Zi = fj must be a singular point for i = 1,2 by Pringsheim's theorem (see, e.g., 
Theorem 17.13 in Volume 1 of Markushevich ifTTl ). In addition to this point, we 
need to find all singular points on \z\ = f, to get the tail asymptotics as we will see. 
As expected from Lemma [6] z = — fj may be another singular point, which occurs 
only when (v-a) does not hold. 

We focus on these singular points instead of searching singular points on \z\ = fj, 
and show that there is no other singular point on the circle through analytic behavior 
of 9i(z). Since results are symmetric for <pi (z) and (piiz), we only consider <p\ (z) in 
this section. 

For this, we use the stationary equation < tl3J >, which is valid on §!. Plugging 

(zi,z 2 ) = (z,? 2 (z)) into Q yields, for |z| e $' mm \fi), 

, (fe(z,yz))-l)fc(C 2 (z)) (Tb(z,C 2 (z))-l)v(0) 

<Pi (z ) = f - -. r . ~ 2 + -I , j . — • (29) 
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Fig. 3 The shaded area is y~ (- ' max) ) nSf+(w ( , 1 ' max) ) 




In the light of this equation, the dominant singularity of <pi(z) is caused by C 2 ( z )' 



(30) 



In addition to %(a,b), we will use the following sets for considering analytic re- 
gions (see Figure[3]l. 

Cs(u) = {z G C;m — 8 < \z\ < u + 8,z ^ u}, u,8 > 0, 



fs(u)=%(u-u ( i' miD) )nc s (u), 



u < —u 



(l.min) 



Remark 5. One may wonder whether ( fT8j ) in Lemma [5] is sufficient for verifying 
analyticity of <pi(z) in Sfj"(Kj max ^) when fi = Hj ™ 1 . This will turn out to have 
no problem because of (|29|). 



In what follows, we first consider the case when (v-a) holds, then consider the 
other case. 



3.3.1 Singularity for the non -arithmetic case 



Assume the non-arithmetic condition (v-a). We consider the analytic behavior of 
<Pi (z) around the singular point z = t\ . This behavior will show that there is no 
other singular point on \z\ = fi. We separately consider the three causes which are 
discussed above. 



([3^) The solution of ( 30 »: This equation has six solutions at most because it can be 



written as a polynomial equation with order six. z = 1 , u[ { are clearly the solutions. 
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Because <f>i(z) of (29i is analytic for \z\ < fi, (30l can not have solution such that 
\z\ < f\ except for the points where the numerator of the right hand side of ( 29 1 
vanishes. This must be finitely many because the numerator vanishes otherwise by 
the uniqueness of analytic extension. On the other hand, ( 30 1 has no solution on the 
circle \z\ = u^'^ except for z — u[ 1 '^ by Lemmajs] 

Thus, the compactness of the circle implies that, if Z\ = < u^' m , then 



( 30 1 has no solution on Cg(i4 ) for some 8 > 0. Hence, we have the following 
fact from d29b. 



Lemma 9. Assume that f\ = u\ ' r ' < u\ ,max) and ^(^(z)) is analytic at \z\ = u\ . 

fir) (I r) 

Then, (p\ (z) has a simple pole at z — u\ , and analytic on C§ (u\ ). 

Remark 6. For categories I and HI, the analytic condition on <P2(C 2 ( Z )) m l emma 
is always satisfied because Lemma|8jand the category condition, C 2 (^i ) < %2> imply-. 

for |z| = u[ l ' T \ 

\m 2 (z))\ < ui^m < m 2 (\z\)) = < o-. 

If fi = wj 1 '^ = u[ l ' max \ then the analytic behavior of <pi (z) around z = Mj is a 
bit complicated because C 2 (z) is also singular there. We will consider this case in 
Sectiong] 

The singularity of C 2 (z): ^ v Lemma |8j this function is analytic on 
^(u^Kuf'^) and singular at z = Kj max \ which is a branch point. 
(|3j:) The singularity of <P2(C 2 ( Z )) : This function is singular at z = fi if C 2 (^i) = 
Otherwise, it is singular at z — u^' m because £ 2 (z) is singular there. Further- 
more, we may simultaneously have £ 7 (ti) = f 2 and fi = u[' max . Thus, we need to 
consider these three cases: fi = u^ ,mw ^ for categories I and HI, and fi < Mj l max ' 

or fi = Kj 1,ma3 ^ for category II. For this, we will use the following fact, which is 
essentially the same as Lemma 4.2 of |[2"T1 . 

Lemma 10. £, (e e ) is a concave function of 9 e [e [ 2 Mn) , 2 (2 ' max) ], ti (u ( 2 l ' max) ) = 0, 
Ci'(4 Lmax) )<0,and 

hm — T ~ z = (31) 

(M (l-x )_ z) . , „ 

Proof. The first part is immediate from the facts that r + is a convex set and i/j 1,max ' = 

(l,m a x) . -p , . (l.max) . (2,max) 

<? D i . By Taylor expansion of ^ 1 (z2j at zi = u\ < u \ 
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~p , \ (l.max) , 1"*^"/ (l,max)w (l.max) \ 9 , /i fl,max)|9\ 

fl(Z2)=«l +^Cl(«2 )(Z2-«2 ) +0(|Z2"«2 I)' 



Letting zi = C 2 (z) m ^ s equation yields ( jjTj ) since C 1 ( C 2 ( z ) ) = z f° r z t0 ^ e su ffi~ 
ciently close to «j I □ 

Another useful asymptotic is: 
Lemma 11. If h! 1, ^ = mS 1 ' , then for any 5 > 0, 



/ (l.max) ,1 J >-"/ (l,max)x 
nm , - } 1. , X s = 777 TT--, ■ (32) 



z— ►« 



Proof. By the condition j^ 1 '™*) = Wj we have 

i-rifcC 2 W)-ri(u (Lmax) )-7ifcC 2 W) 

= (4 1,mo) -i a w)rfi , (»f' ) ) 

+C 2 W(^ ) (« ( 1 1 ' r) )-^ ) W)+^ ) (" ( 1 1 ' r) )-^ ) W. 

Hence, by dividing both sides by (u[ 1,max ^ — z) 1 , Lemma 10 yields |3l| ) because 
p^l (z) and pty (z) are analytic except for z — 0. □ 

We now consider the three cases separately. 

C 2 (^i) < ^2, equivalently, categories I or HI , and f] = itj 1 '™"^: In this 
case, ^(z) is analytic for z € C l 5(M 2 1 ' max ' > ) for some 5 > because « 2 1 ' max - > = 
£ 2 („( 1 - max )) = C 2 (t0 < f 2 . Hence, by Taylor expansion, we have, for \z\ < £2, 



~ / \ ~ / l.max K , -// l.max w ll.maxK , ( \ i.maxlix 

<?>2(Z) = <p 2 04 + %(«2 )(Z- M 2 )+°(|z-"2 I)- ( 33 ) 

(z)) around z= u\ 



Thus the analytic behavior of ^i(^Az)) around z = «|' is determined by that 



of C 2 (z) — « 2 1,max '. Since i^ 1 '™"^ = ^ 2 («[ 1,max ') < f 2 by the conditions of 
Lemma [10] yields 

~ 1 1 (l.max) \ 

- rr 1 w ~i (l.maxK V £<p2\ u 2 > 1 (l.max) >J 
<P2(C 2 (Z)) = <P2(« 2 0- / = , =r("l Z )~ 

/ ~7" ( l.max) s 

+ (|z-4 1 ' max) |z). (34) 

Thus, ^2(C 2 (z)) has a branch point of order 2 at z = fi = Hj 1,max \ and is analytic on 
#+( M ( 1 1,max) ) for some <5>0. 
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(3 



2) £ (f i ) = T2 and X\ < u\' . This is only for category II. Hence, T2 = uk < 



i<2 maX \ an d therefore <p 2 -version of Lemma0is available. Thus, (friz) has a simple 

(2 r) (2 r) 

pole at z = u\ ' ' . Here, u\ is the solution of the equation: 

7 2 (C 1 W,z) = l (35) 
is crucial. Furthermore, £.(z) is analytic at z = fj. Hence, (feiC (z)) has a simple 



12 

pole at z = fi, and is analytic on Cs(u^ ,max ^) for some 5 > 0. 

( 3^-3) C 2 (fi) = fe, and fi = Mj 1,max ': This is also only for category II. This case 

is similar to (3j;-2) except that C 2 ( z ) nas a branch point at z = t\ = u[ 1,max K Since 
<Pi{z) has a simple pole at z = %%, we have, by Lemma 10 



k / v t \ \ r ( 1 .max) \ — l 
<P2(£ 2 (z)) ~ ("i -Z) 2 , 

and (f>2(C 2 ( z )) * s analytic on (i/ 1 1 ' max -') for some 5 > 0. 



3.3.2 Singularity for the arithmetic case 



We next consider the case that (v-a) does not holds. That is, the Markov additive 
process for the interior is arithmetic. In this case, the singularity of <pi (z) at z = f i 
occurs similarly to in Section 3.3.1 In addition to this singular point, we may have 
another singular point — fi as can be seen in Table[T] For this, we separately consider 
two sub-cases: 



(Bl) either (v-b) or w>}> = holds. (B2) neither (v-b) nor mV> = holds. 



(i) 



In some cases, we need further classification: 

,(2) 



(CI) either (v-c) or m = holds. (C2) neither (v-b) nor m 



(2) 



holds. 



Consider (Bl). From Table fl] the solutions of |l6jl are z = ±Mj 1,max ', and the 

J ~~[ i 1 r) 

solution of ( 19 1 is z = u\ . There is no other solution. We consider cases similar to 
(f3}i), (j3>), {3^2), <[3}>1) and ([3}>3) of Section |3~3TT] 

(|3^') The solution of pO"} : This case is exactly the same as in Section 3.3.1 be- 
cause z = — Mj l r ' is not the solution of ( |l9| . Hence, Lemma|9]also holds true. 



([3J3' ) The singularity of £ 2 (z) at |z| = Kj m . It is singular at z — ±u l 

(|3|;') The singularity of <P2(C 2 (z)) a * kl = ^i- F° r z = T\, the story is the same as 
in Section 3.3.1 Hence, we only consider the case that z = — f 1 . From ( 15 1 and 
the condition that (v-a) does not hold, we have 



(l.max) 



2(p 



1 _ POO - y , - s 

— Z2VTl = -C 2 T! , 

-ll+PllTf) 



(36) 
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Hence, |C 2 (-Ti)| = C 2 (Ti) > 0, and 

Since C, (C 9 (^i)) < <Pi(C, (C 9 W)) i s analytic around z=—x\. Furthermore, 



Lemma 10 and ( 34 1 are still valid if we replace uf' m2y ^ by — i<j 1:max ' for i = 1,2. 
However, this z = — f] can not be the solution of (f30|> because of (Bl). Thus, 



we have to partially change the arguments in Section 3.3.1 



l[3j;'-l) C 2 (^0 < ^2 and fi = «j This is only for categories I and HI, and 
q>2(£ 2 (z)) has a branch point of order 2 at z = — u[ l ' max \ and is analytic on 
< ^7{— u^[' ms ^) n ^^(mi 1 '" 13 "') for some 8 > because it has also a branch 

• ^ ^ (l.max) 

point at z = u\ 

|3};'-2) C 2 (^i) = ^2 an( J < »i This is only for category II. Since £ 2 (z) is 
analytic at z = fi, <P2(C 2 W) i s analytic at z = — fi if (CI) holds. Otherwise, 
if (C2) holds, it has a simple pole at z = — fi because C 2 (~ ^l) = ~C 2 (^i) 
is the solution of (|33" 



'-3) C 9 (^i) = ^2 and fi = M| lmax ': This is only for category II, and the 
situation is similar to (|3j;'-2) except that the singularity is caused by C, (z) 
at z = — 1\. To verify this fact, we rework on <P2(C 2 W)- Similarly to (29 1, 
we have, for |z| g (i< 2 ' mm ',f2), 



9z(z) 



(m^-m^iz)) ()b(c 2 ( Z ), Z )-i)v(o) 



Substituting £ (z) into z of this equation, we have 



fc(C,(z)) = 



i-MCj^W),^)) 
(yo(C 2 (C 2 W),C 2 W)-i)v(o) 
i-fcCC^z)),^)) ' 



(37) 



By the assumptions of ([3j;-3), if (C2) holds, then <f>2 (z) has a simple pole at 
z = — Vi, and therefore (^(^(z)) ~ (— «j 1,max ' — z)~z around z = — Mj 1,max ' 
by LemmajToj Otherwise, if (CI) holds, we need to consider <Pi (Cj (C 2 C0)) 
in ( Jy/j i due to the singularity of £ 2 (z) at z = — fi = — M| 1,max ', where 
(piiC^z)) is analytic atz = — « 2 1 ' max ' ) = — £ 2 («[ 1 ' max ') because 

\z l (c 2 (-ti))\=\z 1 (t 2 m\<tt- 
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Hence, ft (C 1 (C 2 (z))) - ft (-£ (w^)) ~ (-" ( /' max) -z) * ■ On the other 
hand, 7i(Cj(C 2 (z))>C 2 (z)) _ 1 ~ (-«^ max ' ) -z) J because (v-a) does not 
hold. Combining these asymptotics in (37 1, we have (piiC 2 ( z )) ^ ^(^^2) ~ 

(— m^ 1 '™"' — z) 2 around z = — z^ 1 '™"^ by Lemma [To] 
We next consider (B2). From Table |1| the solutions of il6\ are z = ±Wj max \ 



and the solutions of jl9| ) are z = ±m[ 1i; . In this case, the arguments for z = — ft 
is completely parallel to those for z = f\ except for the cases ([3};'-2) and (|3j;'-3). 
The latter two cases are also parallel if (C2) holds. Otherwise, 92 (z) is analytic at 

Z = — T2. 



3.4 Asymptotic inversion formula 

From these singularities, we derive exact tail asymptotics of the stationary distribu- 
tion. For this, we use Tauberian type theorem for generating functions. 

Lemma 12 (Theorem VI.5 of [7|). Let / be a generating function of a sequence 
of real numbers {p(n);n — 0,1,...}. If /(z) is singular at finitely many points 
a\,a2, ■ ■ ■ ,a m on the circle |z| = p for some p > and positive integer m, and ana- 
lytic on the set 

Ai = {zE C;|z| < r u z ^ a;, |arg(z-a,-)| > i = l,2,...,m, 
for some 0), and r ( - such that p < r, and < CO, < f and if 

lim (ai-z) Ki /(z)=^, i=l,2,...,m, (38) 

for JQ ^ {0, — 1 , —2, . . .} and some constant b, E R, then 

E iFR «rj p(»)=6, (39) 

for some real number b, where F(z) is the gamma function for complex number z 
(see Sec 52 of Volume II of ifTTl ). 

Recall that the asymptotic notation "~" introduced in Section[T] With this nota- 
tion, ( |39l ) can be written as 

where T(i) = 0r andT(-i) = -2-s/n. 
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We will apply Lemma[l2|in the following cases: For m= l,a\ = Mj 1 '^ and K\ — 

1 ^ (l.max) ! , 1 ,-. , (l,r) * , ~ , M.max) 

\,2,ai=u\ and K\ = ±k. For m = 2, a\ = ±u\ and fCi = 1,2, ai = ±u\ 
and )fi - 1 



2- 



4 Exact tail asymptotics for the non-arithmetic case 

Throughout this section, we assume the non-arithmetic condition (v-a). We first 
derive exact asymptotics for the stationary probabilities v(n,0) and v(0,n) on the 
boundary faces. Because of symmetry, we are only concerned with v(n,0). 



4.1 The boundary probabilities for the non-arithmetic case 

(if) (2 r) (i r) (2 r) 
We separately consider the two cases that u\ ' <u 2 ' and u 2 ' ~>u\' , which 

correspond with categories I (or III) and n, respectively. In this subsection, we prove 

the following two theorems. 

Theorem 2. Under the conditions (i)-(iv) and (v-a), for categories I and M , fi = 
(ID 

u\ ' , and P(L\ = n,L,2 = 0) has the following exact asymptotic h\ (n). 

/M(»>=< n-U-\ = "! Lmax) '= " ( r r) , (40) 

-I — n (IT) (l,max) , (l,r) 

By symmetry, the corresponding results are also obtained for P(L\ — Q,L2 — n) for 
categories I and II. 

(2 r) 

Theorem 3. Under the conditions (i)-(iv) and (v-a), for category n, ?2 = u 2 ' , and 
P(L\ = n,L,2= 0) has the following exact asymptotic h\ («). 



(l,max) (l,r) 





< Mj 


, or 


?1 


(IT) 


= u 


?1 


(IT) 


7^ U 


?1 


(IT) 
— u\ 1 


= u 



*i(») = < rinTl^r" 1 ' (41) 



= u\ 

(l.max) 

H , n — mj 7- wj 

-i=-n * (l.max) , (l,r) 



By symmetry, the corresponding results are also obtained for P{L\ — } L2 — n) for 
categories HI. 



Remark 7. Theorems [2] and [3] are exactly corresponds with Theorem 6.1 of [4| (see 
also Theorems 2.1 and 2.3 of [3]). This is not surprising because of the similarity of 
the stationary equations although moment generating functions are used in SHI- 



24 M. Kobayashi and M. Miyazawa 

Remark 8. These theorems fill missing cases for the exact asymptotics of Theorem 
4.2 of [ 18 1. Furthermore, they correct two errors there. Both of them are for category 
II. The exact asymptotic is geometric for fi = Kj r ' = m^ 1 '" 13 "' = Kj , and not 

geometric for X\ = u^' 1 ' / u^ ,msa ^ (see Theorem 3). However, in Theorem 4.2 
of [18], they are not geometric (see (43d3) there) and geometric (see (4c) there), 
respectively. Thus, these should be corrected. 

Proof of Theorem |2j^ We assume that either in category I or IE occurs. This 

is equivalent to w 2 ' ^ < u^'^, and Ti = u\ ,r \ Furthermore, we always have 

f i r) i i r) i l r) 

C 2 ( M i ) = u 2 < an d therefore <P2(C 2 W) * s analytic at z = u\ ' . We con- 
sider three cases separately. 

Qi) u[ ir) < u[ h ™ x] : This case implies that « ( /_ r) < M [ 1,max) and fi(u^' m ^) > 1, 
and therefore . Hence, by LemmaM <p\ of d29k satisfies the conditions 



12 



under the setting (38 1 with a\ = u\ , K = \ . Thus, letting 

b = (fc(n^) - 1)92(4' J) ) + (ftCggg) - 1 ) y (°) 
s7i( M 'C 2 («))l l(=H (i. r ) 



(1 r) 

which must be positive by (39 1 and the fact that <pi(z) is singular at z = Mj , we 
have 

limf!'P(Li =n,L 2 = 0)=b. 

(J^Jd) Mj 1,r ' = Mj 1,max \ Mj 1 ' 1 ' = Mj l max '': In this case, category HI is impossible, and 

fi (u( lmax )) = 1. On the other hand, <pz(z) is analytic atz= ^ 2 {u^' m3X ^) < T 2 because 
of Category I. Hence, we can use the Taylor expansion ( [33] ^ and therefore ( 29 1, ( 34 1 
and Lemma[TT|yield, for some 5 > 0, 

lim (u[ 1 >™ x) -z)hi(z)=b, (42) 

,.sr_i_ , ( 1 ,max) ( 1 ,max) 

where 



((r 2 (u( 



1)^2(4™) + (]b(ul 1 «)-l)v(0) 



V2P1W) 



'1 



Hence, <pi satisfies the conditions of Lemma 12 under the setting ( |38] > with a 
Mj max ' and JCi = j, and therefore we have 

]imnh?P(Li =n,L 2 =0) = -4=. 

where the positivity of b is checked similarly to case (|4^) (see also case (j4f;) below). 
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([4]:) Hj 1 '^ = Mj 1,max \ u[ 1 '^ ^ «[ 1,max ': In this case, category HI is also impossible, 
and y 1 (u( 1,max )) ^ 1. Thus, we consider the setting (38 1 with K\ = — \. From (29 1, 
we have 

<^)-<Pi(" ( i 1,maX) ) 

(f 2 (z, C 2 W) - l)fc(Cj(z)) + (ft(z, C 2 (*)) - l)v(O) . (1 max , 

= — " <Pl( " 1 } 

= (fefe C 2 W) - iKfeggW) - fe(4 1,max) )) 

(fe(z,C 2 (z))-fe(u( 1 - mm )))fe(4 1,m " ) ) (yofeC 2 W)-yo(u (1 - max )))v(Q) 

+ l-ri(z,C 2 W) + \-%{z,^{z)) 

(fe(u (l,max) ) _ 1) ^ (M (l,max) ) (y ( U ( LmaX >) - l)v(O) 



(i-M^C 2 M))(i-yi(u(«)) (i-yifei 2 (z)))(i-yi(u( 1 ' m -))) 

x(7ifcC 2 W)-7i(u (1 ' max) )). (43) 
We recall ((34]) that 

~ 1 { (l.maxK 

~ / ^ - 1 (l,max)s / (l.max) \ \ V 1< ?2\ U 2 ) , i\ (l.max) ,J, 
<P2(C 2 W)-<P2(« 2 ') = -(»l -Z) 2 ^== , + °(K ~Z| 2 )- 

From pTj ), we have 

ftfe^z)) - 7o(u (Lmax) ) = (C 2 (z) - C 2 (4 1 ' max) ))/>i°i ) W 



(0)/ (l.max) n 

>/2 Ptl ( Hl ) (M ( 1 , m ax)_ z) . +o(|M ( 1 , m ax)_ z| ., 

Ci(4 1,max) ) 



Similarly, 

7i(z,C 2 «)-7i(u (1 ' max) ) 

at (1)/ fl,max)\ 

^* l(Ml \u^-z)^+o(\u^-z\h, 



72(z,C 9 «)-72(u (Lmax) ) 



(2) / (l,max)\ 
I ' 1 ■ / (l,max)\2 



-r 1 '(4 i,max) ) 



/ ( 1 ,max J \ z 1 

\ U 2 I / 1 (l.max) , /1 (l.max) i'\ 

— —{u\ ' -zY~ +o{\u\ 
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With the following notation, 



c\ = 



V2 



-c'i'(4 1 ' raax) ) 



d k 



dv 



(l.max) 
1 



(l.max) 



43b yields, as z 



M (l,max) sat j s fyj n g t jj at z g &g(u 



( 1 ,max) 
1 



for some 5 > 0, 



- ^(« ( /' max) ) = -ciCk?'— J ((r 2 (u (1 ' max) ) l)^(4 1,max) ) 



+rf 2 9 2 (4 1 - max) ) + *v(0)+./ 1 ^ 1 (« ( 1 Umax) ))+ o (|« ( 1 1 ' raax) -z|^). 



(44) 



Let 



Then, taking mi which is sufficiently close itj m from below in (|44|), we can see 

that this b must be negative because <j>i (u\) is strictly increasing in u\ e [0, 4 

1 

2' 



(l,max) 



Thus, <|3S b holds for the setting of <|3 8 b with K\ = —\, and therefore (|39| leads to 



lim«3ffP(Li = n,L 2 =0) 



20r 



>0. 



Thus, we have obtained all the cases of (40 1, and the proof is completed. □ 

Proof of Theorem^ Assume category II . In this case, T 2 = £ 2 (?i), and <p 2 (z) has 

a simple pole at z = f 2 because of category II (see (|3]>2)). We need to consider the 

following cases. 

( i r) 

(4i'): fi < u\ ' : In this case, <p 2 (C 2 (z)) has a simple pole at z = t\- Since <Pi(z) 

has no other singularity on |z| = fi.it has a simple pole atz = f\. 

Up'): fi =H ( 1 1 ' r) :This case is further partitioned into the following subcases: 

(|4]3'-1) Mj r ' 7^ Mj 1,max ': In this case, fi = 4 < Mj 1,max \ and therefore it is easy 
to see from ( 29 1 that <pi (z) has a double pole at z = f i . Hence, we can apply the 



setting <|3 8 b with a\ = 4 an d K = 2 



d4b'-2) u 



M 



(l.max) 
1 



7^ M 



( 30 1 does not hold, and therefore |3l| and the 
fact that fyziz) has a simple pole at z = f 2 yield the same asymptotic as ( j42) > but 
with a different b. Hence, we apply (38 i with Ki = j. 

(|4j?' -3) 4 1 ' r ' = Mj 1,max ^ = 4 1 • m this case, we note the following facts . 
(|4]y-3-l) <p 2 (z) has a simple pole at z = f 2 , and therefore Lemma 10 yield 



fc(C,(z))~ (« 



(l.max) 
1 



-z)" 
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gJ>'-3-2) ByLemma[n] 1 - fi (z,C 2 (z) )-(«', 
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l.max) 



■Z 2- 



lr, , , 

«i = Mj and JC = 1. 
Thus, similar to Theorem|2j we can obtain (41 1, which completes the proof. □ 



Hence, (29 1 yields (p\{z) ~ (mj max; — z) 1 , and therefore we apply (38 i with 



4.2 TTie marginal distributions for the non-arithmetic case 



We consider the asymptotics of P((c,L) > jc) as jc — > °° for c = (1,0), (0,1), 
(1,1). For them, we use the generating functions p+(z, 1), <p + (l,z) and (p + (z,z). 
For simplicity, we denote them by l//io(z), Vbl(z), Vn(z), respectively. We note 
that generating functions are not useful for the other direction c because we can 
not appropriately invert them. For general c > 0, we should use moment generating 
functions instead of generating functions. However, in this case, we need finer an- 
alytic properties to apply asymptotic inversion (e.g., see Appendix C of 0). Thus, 
we leave it for future study. 
From i 



<p(z) 



13 1 and (28 1, we have, for z g C 2 satisfying (\zi |, \zi\) S 
7i(z) 



1 + 



1 



!-?+(*) 



ft(zi) + ^^ fefa) + M^^ v(0). (45) 



!-?+(* 



Hence, the asymptotics of P((c,L) > x) can be obtained for c = (1,0), (0, 1), (1, 1) 
by the analytic behavior of y/io(z), ^oi (z), V^ii (z), respectively, around the singular 
points on the circles with radiuses p c , where 

P(i,o) = sup{n>0;(w,l) e p ( o,i) =sup{«>0;(l,n) € #}, 
P(i,i) = sup{w > 0; (w,w) € ^}. 

Since y/io(z) and y/bi(z) are symmetric, we only consider yio(z) and Vn(z). 
From ( |45| ), we have 



Vio(z) = 1 



Yn(z) 



l-f+(z,l) 



7i(z,z)-l 
l-7+(z,z) 



+ 



l-7+(z,l) 
jb(z,l)-y+(z,l) 

l-7+(z,l) 



v(0), 



- , x . 72(z,z)-7+(z,z) - , , 

<Pl (Z) + : r 92 (z) 



+ 



l-7+(z>z) 
?o(z,z)-?+(z,z) 
l-7+(z,z) 



v(0). 



(46) 



(47) 



We first consider the tail asymptotics for c = (1,0) under the non-arithmetic con- 
dition (v-a). From (46 1, the singularity of yfu(z) on the circle \z\ = P(i,o) occurs by 
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either that of (f>\ (z) or the solution of the following equation: 

f+(z,l) = l. (48) 

Since this equation is quadratic and the domain <%! contains vectors x> 1= (1,1), 
the equation (48 1 has a unique real solution greater than 1. We denote it by o+. We 
then have the following asymptotics (see also Figure[4]). 

Theorem 4. Under the conditions (i)-(iv) and (v-a), let h\{n) be the exact asymp- 
totic function given in Theorems [5] and [3] then P{L\ > n) has the following exact 
asymptotic g\ («) as n — > oo, 

(a) If C 2 (" ( i LF) ) < 1. thengi(n) = a+" . 

(b) If C 2 («i 1,r) ) > 1 and i 2 (« ( i 1,r) ) # !> thengi(n) =Ai(n). 

(c) If C 2 (" ( i Lr) ) > 1 = C 2 (" ( r r) ), then^n) = ff". 

(d) If C 2 (" ( /' r) ) = 1 = C 2 (" ( r r) X then ft (n) = Tf". 

(e) If C 2 (" ( i Lr) ) = 1 > C 2 (" ( r r) ), then ft (n) - nff". 

Remark 9. The corresponding but less complete results are obtained using moment 
generating functions in Corollary 4.3 of [ 18 1. 

Before proving this theorem, we present asymptotics for the marginal distribution 
in the diagonal direction. Let Od be the real solution of 

7+(m,m) = 1, u > 1, 

which can be shown to be unique (see Figure [6j. Because of symmetry, we assume 
without loss of generality that fi < f 2 . See Figure[5]for the location of this point. 

Theorem 5. Under the conditions (i)-(iv), (v-a) and fi < f 2 , let h\(n) be the ex- 
act asymptotic function given in Theorems [5] and [3] then P(L\ +L 2 > n) has the 
following exact asymptotic g+{n) as n — > oo. 

(a) If a d < fi, then g+(n) = a d ". 

(b) If cr d > fi, theng + (n) =h\(n). 

(c) If a d = fi 7^ M ( 1 1 ' max) , then g + (n) = na d ". 

(d) If C7 d = fi = M ( j 1,max) = f 2 , then g + (n) = na d " 

(e) If a d = fi = M^' max) ^ f 2 , then g+(n) = a"". 

In what follows, we prove Theorem[4] The proof of Theorem[5]is similar, so we 
only shortly outline it. 
Proof of Theorem^ Let 

S(z) = {Hz, 1) - ?+(z, l))fc(l) + (%(z, 1) - f+fe l))v(0), 



then |46) can be written as 
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Since 1) > !,%(«, 1) > 1 for m > and 



ft 0,1) 



<0, 



U=0\ 



du 



f+(u,l) 



-f+(z,iy 



>0 ifC(cr + )=0, 



(49) 



W=<7-f 



where Ci is a positive number satisfying that f\ ((7i, 1) = 1, % (o+) > and o+ = G\ 
implies that the prefactor of <p\ (z) is positive at z — <J+ if C 2 (°+) = 0- Having these 
observations in mind, we prove each cases. 

(a) Assume that C2( M j_ ) < 1. This occurs if and only if a + = pio < t\ (see the 
left picture of Figure H|. In this case, ywiz) must be singular at z — <J+ because it 
is one the boundary of the convergence domain @. Hence, it has a simple pole at 

z = <7+, and therefore we have the exact geometric asymptotic. 

— ft r) l i r) 

(b) Assume that £ 2 ( M i ) > 1 and C 2 ( M i ) !• This case occurs if and only 

if a + 7^ pio = fi (see the right picture of Figure [4J. In this case, fi(fi, 1) 7^ 1, 
7+(fi, 1) 7^ 1 and fi(fi, 1) — 1 has the same sign as 1 — f+(ii, 1). Hence, the pref- 
actor of (pi(z) is analytic at z — f\, and the singularity of ty/io(z) is determined by 
q>i (z). Thus, we have the same asymptotics as in Theorems [2] and [3] 
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(c) Assume that £ 2 ( M i ) > 1 and C 9 ( M i ' ) = 1 ( see the left figure of Figure 5 1. 
In this case, y+(fi,l) = fi(fi, 1) = 1 and category II is impossible, and therefore, 
from (49 1 and Theorem|2] we have the exact geometric asymptotic. 

(d) Assume that £ 2 ( M i ) = 1 = C 2 ( M | 1 F ') ( see tne right figure of Figure In 

and therefore f+(zi, 1) = 1. We need to consider two 

thenyi(fi,l) = 1 



, - - (l.max 

this case, i\ = (7+ = u\ 



subcases, = 



(l.max) 
1 



and u 



1 



^ W j ma . If Mi 



U 



(l,max) 
'l 



and (pi (z) ~ (fi — z) 2 by Theorem]^] Thus, we have 1/10(2) ~ (^i — z) 1 due to the 
second term of ( 49 1. Otherwise, if u\ ^ 7^ u[ l ,max ' , then fi (fi , 1) 7^ 1 implies that the 



prefactor of <pi (z) in (46 1 has a single pole at z = fi and that (p\ (z) — <pi («[ 1 ' max ' ) ) 

(fi — z) 1 ■ Again from (49 1, we have y/io(-z) 
geometric asymptotic in both cases. 



(fi — z) . Thus, we have the exact 



(e) Assume that £ 2 ( M i 



) = 1 > C 2 ( M i ' )• In this case, %\=a + = u x 



( 1 ,max) 



and we must have category I or HI. Since f+(fi, 1) = 1, fi(fi, 1) > 1 and <pi(z) has 
a single pole at z = Ti, 1/10(2) in (46i has a double pole at z = Ti. This yields the 
desired asymptotic. □ 



Fig. 6 Left: cr d < f 1 , Right: <7 d > f 1 
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The proof of Theorem[5]is more or less similar to that of P{L\ > n). From Figures 
[6]and[7] we can see how the dominant singular point is located. Since its derivation 
is routine, we omit detailed proof. 



5 Exact tail asymptotics for the arithmetic case 



Throughout this section, we assume that (v-a) does not hold. As in Section 3.3.2 we 
separately consider two cases: (Bl) either (v-b) or = holds, and (B2) neither 
(v-b) nor = holds, according to Table 1 In some cases, we need: (CI) either 

(2) (2) 

(v-c) or m\ = holds, and (C2) neither (v-c) nor m\ = holds. 



5.1 The boundary probabilities for the arithmetic case with (Bl) 

In this case, we have the following asymptotics. 

Theorem 6. Under the conditions (i)-(iv) and (Bl), if (v-a) dose not hod, then for 
categories I and IH , fi = u\' , and P{L\ = n,L,2 — 0) has the following exact 
asymptotic hi(n). For some constant b E [—1, 1], 

l'2(n)= { n-U-\ «j 1,n =« ( /' maX Wi 1,r) > (50) 

n-\ (l + ft(-l)»)ff», M ( /' r) = u[ l ' max) £ u[ hl) . 

By symmetry, the corresponding results are also obtained for P(Li —0,L2 — n) for 
categories I and II. 

Theorem 7. Under the conditions (i)-(iv) and (Bl), if (v-a) dose not hod, then, for 
category II, T\ = £ (T2), ^2 — u\ ' , and P{L\ = n,L2 = 0) has the following exact 
asymptotic hi(n). For some constant b G [—1, 1], 

ff ", fi < u^' or 

(IT) (l,max) (Is) 
Tl = Mj = Mj — U\ , 

— n - (l.-T) / (l,max) 

«Tj , Tl = U\ ' f: u\ 

n-kp, f l=u [^= u f^^ u ^\ (51) 

and (CI) holds. 

and (C2) holds. 

By symmetry, the corresponding results are also obtained for P(L\ =Q,L2 = n) for 
categories HI. 



hz{n) = < 
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Remark 10. As we will see in the proofs of these theorems, the asymptotics can be 
refined for those with the same geometric decay term f j~ ". There is no difficulty to 
find them, but they are cumbersome because we need further cases. Thus, we omit 
their details. 

Remark 11. One may wonder whether b = ±1 can occur in Theorems [6] and [7] If it 
is the case, then the tail asymptotics are purely periodic. Closely look at the coeffi- 
cients of the asymptotic expansion of the terms in ( |37l >, this unlikely occurs because 
\pz(— C 2 (^i))l < 92(C 2 (^i))- Thus, we conjecture that \b\ < 1 is always the case. 

By Table [TJ <pi(z) may be singular at z = — fi on |z| = fi. On the other hand, 
<Pi (z) has the same singularity at z = 1 1 as in the non-arithmetic case, so we can 

only focus on the singularity at z = — fi. We note that z = —u\ ' can not be the 
solution of ( fL9| ) under the assumptions of Theorems [6] and [7] Having these in mind, 
we give proofs. 

Proof of Theorem^ We consid er the singularity of <pi (z) at z = — X\ by (29 1 using 



the arguments in Sections 3.3.2| and|4| Note that fi = u[ 1,F ^ because the category is 



either I or HI. We need to consider the following three cases. 

@i) uf T) ± M ( ! Lmax) : This case is equivalent to u^ ,r ' < u[ 1 ' mRX \ and it follows 
from (29 1 that <pi(z) is analytic at z = —u^ . Hence, there is no singularity 



contribution by z = — Uj r '. 

m x (l.r) (l,max) fl,rl (l.max) T ^, . (l.max) . . 

(|5b) u\ =u\ ,u\ =u\ : In this case, as z — > — Mj in such a way 
that z G &+(-u[ UmRx) ) for some 8 > 0, 



, <* , (l.max) \\ / (l.max) \\ 

<P2(C 2 (z))-<P2(C 2 (-«i -z) 2 

but 1 — 7i (z, C 2 ( z )) does not van ish at z = — Mj 1,max ', and therefore 

~ / \ ~ / fl,max)\ / (l.max) \ ' 

This yields the asymptotic function n _ 2fj~", but this function is dominated by 
the slower asymptotic function n~3f j~ " due to the singularity at z = u\ l,mRX \ 
|5};) Mj = Mj 1,max ', ^ Mj max ^: In this case, the solution of (19i has no 



2 



essential role, so <pi(z) has the same analytic behavior at z = — mj^ 3 "' as at 
z = «j 1,max ' in ([4J;) in the proof of Theorem^ 

Thus, combining with the asymptotics in Theorem[2] we complete the proof. □ 

Proof of Theorem^ Because of category II, T2 = £ 2 (fi), and therefore £ 2 (— Ti) = 
— £ 2 (fi) = — T2 by the assumption that (v-a) does not hold. We consider the singu- 
larity at z = — f 1 for the following cases with this in mind. 
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(|5ji') X\ <M ( /' r) : This case is included in (3|;'-2). Hence, if (CI) holds, ^(^(z)) 

(1 

and therefore (p\ (z) are analytic at z — —u\ ' . Otherwise, if (C2) holds, <P2(C 2 W) 
has a simple pole at z = — w[ . However, in (29 1, <p2(C 2 (z)) nas tne prefactor, 

fir) 

72(2, C 2 (z)) — 1' which vanishes at z = — u\ ' because of (C2). Hence, the pole 
of (/>2(£ 2 (z)) is cancelled, and therefore <pi(z) is analytic at z — — u[ l ' T \ Thus, 

(1 r) 

either case has no contribution by z = —u\ ' . 

(j5j?' ) fi = «j 1,r ' 7^ «j 1,max ': In this case, fi = u[ lj \ If (C2) holds, then (feiz) has 
a simple pole at z — — fi, and therefore as in (|4j3'-3-l), 



~ /r / w / (l.max) \ 



but (|4j3' -3-2) is not the case, and therefore this yields the asymptotic function 
rr^fy" . However, this asymptotic term is again dominated by ff" due to the 

singularity at z = u^' max \ On the other hand, if (CI) holds, then there is no 
singularity contribution by z = — 1\. Hence, we have the same asymptotics as in 
the corresponding case of Theorem[3] 

(|5j:') Ti = Kj 1 ' 1 ^ — u^'" 1 ^: This is the case of |iJ;'-3). As we discussed there, if 

(C2) holds, <P2(C 2 W) ~ ( — «i — z)~? around z = -Mi ■ Because of 
(Bl), there is no other singularity contribution in (29 1, and therefore we also 

have <pi(z) ~ (— i<J 1 ' maX ' ) — z)~z around z = — u[ 1,m \This results the asymp- 
totic n~?f^ n . On the other hand, if (CI) holds, we similarly have <pi(z) — 

~ / fl,max)\ / (l.max) \1 rm • ■ i* ^i A - — I — n ^ 

<j?i (— u\ ) ~ (— u\ — z) 1 ■ This implies the asymptotic n 2 T t . To com- 
bine this with the corresponding asymptotics obtained in Theorem [3] we con- 
sider two subcases. 



u\ ,m = u\ . In this case, the asymptotics caused by z = f 1 is nf ] 
and therefore the asymptotics due to z = —u[ l ' m ^ is ignorable. 
©'-2) u[ l ' max) ? u™: In this case, the asymptotics caused by z = f\ is 
rT 2 f j~ ". Hence, we have two different cases. If (CI) holds, the contribution 
by z = — fi is ignorable. Otherwise, if (C2) holds, then we have additional 
asymptotic term: (— l)"n~ 2 fj~". 

Thus, the proof is completed. □ 



5.2 The boundary probabilities for the arithmetic case with (B2) 



We next consider case (B2). As noted in Section 3.3.2 in this case, <pi (z) is singular 



at z — ±u[ l,v> , and both singular points have essentially the same properties. Thus, 
we have the following theorems. 
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Theorem 8. Under the conditions (i)-(iv) and (B2), if (v-a) dose not hod, then for 
categories I and M , fi = u\' , and P{L\ = n,L,2 — 0) has the following exact 
asymptotic hj(n). For some constants bj E [—1, 1] for ; = 1,2,3, 

r (i+bri-mxr, u [ un ^ u [ l ™ x \ 

h 3 (n) = I n -2(l + A 2 (-l)»)ff», u[ Lr) = u[ L ™ x] = u[ U \ (52) 

ln-?(i+*3(-i)»)^-», «[ Lr) = «[ Lmax V« ( r r) . 

By symmetry, the corresponding results are also obtained for P(L\ —Q,L2 — n) for 
categories I and II. 

Theorem 9. Under the conditions (i)-(iv) and (B2), if (v-a) dose not hod, then, for 
category II, T2 = u^K an d ^(^l = = 0) has the following exact asymptotic 
/i3(n). For some constants bj e [—1,1] for i = 1,2,3, 

(H-Ai(-l)»)ff", fi<«S 1>r) Qr 

(IT) (l.max) (l,r) 
Tl=ll\ =U\ =U\ >, 

n(l+b 2 (-lY)t-", i^u^^u^, ( 
n- l i(l+b 3 (-l)")x- n , z^u^^uf'^^u^. 

By symmetry, the corresponding results are also obtained for P(L\ —Q,L2 — n) for 
categories HI. 




5.3 The marginal distributions for the arithmetic case 

Under the arithmetic condition that (v-a) does not hold, we consider the tail asymp- 
totics of the marginal distributions. Basically, the results are the same as in Theo- 
remsffland Theorem[5]in which Theoremsj2]andj3] should be replaced by Theorems 
[6]and[7]for the case (Bl) and Theorems [8]and [9]for the case (B2). Thus, we omit 
their details. 



6 Application to a network with simultaneous arrivals 

In this section, we apply the asymptotic results to a queueing network with two 
nodes numbered as 1 and 2. Assume that customers simultaneously arrive at both 
nodes from the outside subject to the Poisson process with rate X. For i — 1,2, ser- 
vice times at node i are independent and identically distributed with the exponential 
distribution with mean pL^ 1 . Customers who have finished their services at node 1 
go to node 2 with probability p. Similarly, departing customers from queue 2 go to 
queue 1 with probability q. These routing is independent of everything else. Cus- 
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tomers what are not routed to the other queue leave the network. We refer to this 
queueing model as a two node Jackson network with simultaneously arrival. 

Obviously, this network is stable, that is, has the stationary distribution, if and 
only if 

—, < Mi > —, < M2- (54) 

This fact also can be checked by the stability condition (iv). 

We are interested in how the tail asymptotics of the stationary distribution of 
this network are changed. If p = q = 0, this model is studied in |8]|9). As we will 
see below, this model can be described by a double QBD process, and therefore we 
know solutions to the tail asymptotics problem. However, this does not mean that the 
solutions are analytically tractable. Thus, we will consider what kind of difficulty 
arises in applications of our tail asymptotic results. 

Let L((t) be the number of customers at node i at time t. It is easy to see that 
{(L\(t),L2(t));t £ W.+} is a continuous time Markov chain. Because the transition 
rates of this Markov chain are uniformly bounded, we can construct a discrete time 
Markov chain given by uniformization, which has the same stationary distribution. 
We denote this discrete time Markov chain by {L„ = (Ln,L2i);£ € where it is 
assumed without loss of generality that 

A + jUi +ju 2 = 1. 

Obviously, {L„;£ 6 Z + } is a double QBD process. We denote a random vector 
subject to the stationary distribution of this process by L = (Li,Li) as we did in 
Section[2] 

For applying our asymptotic results, we first compute generating functions. For 

u = (u.\,U2) € M 2 , 

y+(u) = Xu\it2 +flipw[ l U2 + pL2qu\ii2 1 +/Xi(l — p)u~[ l +M2(1 — l) u 2 1 ' (55) 
7i(u) = Xu\U2 + y.\pu\ x 112 + -p)u\ x (56) 
j^(u) = Xu\U2 + Il2qu\ii2 1 +M2(1 -<?)«2 1 (57) 

We next find the extreme point u' 1 ' r - ) = (hj 1 ,1*2 ). This is obtained as the 
solution of the equations: 

7+(u) = fi(u) = l 
Applying ( |55j ) and ( f56*} to the first equation, we have 

U2 = u\q + (l — q). (58) 



Substituting (58 1 into 71 (u) = 1, we have 



Xu\{u\q+ 1 — q) + H\p{u\q + 1 — #) + jUi(l — p) +jU 2 «l =U\. 
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Assume that q > 0. Then, u\ has the following solutions. 

-X±yJX 2 +4A<?jlti (1 - pq) 



Ml = 1, 



2Xq 

l.n 



We are only interested in the solution u\ > 1, which must be u\ ' , that is 



,(l,r) _ -A + y/A 2 + Akq jli(\ - pq) 
2Xq 



u\ L ' T> = — v — ^ , (59) 



We next consider the maximal point u' l max ' of y(u) = 1. This can be obtained 
to solve the equations: 

~ i \ 1 du\ 
7+(u) = 1 ' ^=°- 

These equations are equivalent to 

Ami + \l\pu\ X — p^quxu^ 1 — ^(1 — q)u^ 2 = 0, (60) 



Xu\U2 + \l\pu x l U2 + pl2qu\u 2 1 +jlti (1 — 1 +^2(1 ^ ?)«2 1 = 1- (61) 



Theoretically we know that these equations have two solutions such that u > 0, 
which must be u' l min ' and u' 1 ' max ' ) . We can numerically obtain them, but their an- 
alytic expressions are not easy to get. Furthermore, even if they are obtained, they 
would be analytically intractable. 



Fig. 8 Effect of the arrival rate: A is changed from 1 to 1.2 and 1.5 (thicker curves) while p.i = 5, 
= 4, p = 0.25, q = 0.4 are unchanged 
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To circumvent this difficulty, we propose to draw figures. Nowadays we have ex- 
cellent software such as Mathematica to draw two dimensional figures. Then, we can 
manipulate figures, and may find how modeling parameters change the tail asymp- 
totics. This is essentially the same as numerical computations. However, figures are 
more informative to see how changes occur (see, e.g., Figure[8]l. 

We finally consider a simpler case to find analytically tractable results. Assume 
that q = but p > 0. q = implies that 

u^MprM), 

where p\ = j^. Obviously, p\ must the decay rate of P(L\ > n). This can be also 
verified by Theorem[4] However, it may not be the decay rate of P(L\ > n,L? = 0). 
In fact, we can derive on the curve y+(u) = 1, 



du2 



Hence, u} = Mj ' if and only if 



- Ol ±_hp) 



Mi-pr 1 ) 



M2>Mi+Api. (62) 



Thus, if (62 1 holds, then P{L\ > n,L,2 = 0) has the exactly geometric asymptotic. 



Otherwise, we have, by Theorem|2j 



\\mn- i 2(u ( ^ xn ™ ) )- n P(L x >n,L 2 =0)=b. (63) 

n— s-oo 

We can see that pr 1 < M[ l max ', but Mj l max ' is only numerically obtained by solving 



60) and |6T}. 



7 Concluding remarks 

We derived the exact asymptotics for the stationary distribution applying the analytic 
function method based on the convergence domain. We here discuss which problems 
can be studied by this method and what are needed for further developing it. 

(Technical issue) In the analytic function method, a key ingredient is that the 
function C 2 M ' s analytic and suitably bounded for an appropriate region as 
we have shown in Lemma|8j For this, we use the fact that C 2 ( z ) ^ s the solution 
of a quadratic equation, which is equivalent for the random walk to be skip free 
in the interior of the quadrant. The quadratic equation (or polynomial equation 
in general) is also a key for the alternative approach based on analytic extension 
on Riemann surface. If the random walk is not skip free, it would be harder to 
get a right analytic function. However, the non skip free case is also interest- 
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ing. Thus, it is challenging to overcome this difficulty. One here may need a 
completely different approach. 

{Probabilistic interpretation) We have employed the purely analytic method, and 
gave no stochastic interpretations except a few although the asymptotic results 
are stochastic. However, probabilistic interpretations may be helpful. For ex- 
ample, one may wonder what are probabilistic meaning of the function £ and 



the equation (29 1. We do believe something should be here. If they are well an- 
swered, then we may better explain Lemma [8] and may resolve the technical 
issues discussed above. 

(Modeling extensions) We think the present approach is applicable for a higher 
dimensional model as well as a generalized reflecting random walk proposed 
in |fl9l as long as the skip free assumption is satisfied. One may also consider 
to relax the irreducibility condition on the random walk in the interior of the 
quadrant. However, this is essentially equivalent to reducing the dimension, so 
there should be no difficulty to consider it. Another extension is to modulate the 
double QBD or multidimensional reflecting random walk in general by a back- 
ground Markov chain. The tail asymptotic problem becomes harder, but there 
should be a way to use the present analytic function approach at least for the 
two dimensional case with finitely many background states. Related discussions 
can be found in [ 19 1. 

(Applicability) As we have seen in Section |6| analytic results on the tail asymp- 
totics may not be easy to apply for each specific application because they are 
not analytically tractable. To fill this gap between theory and application, we 
have proposed to use geometric interpretations instead of analytic formulas. 
However, this is currently more or less similar to have numerical tables. We 
here should make clear what we want to do using the tail asymptotics. Once a 
problem is set up, we may consider to solve it using geometric interpretations. 
Probably, there would be a systematic way for this not depending on a specific 
problem. This is also challenging. 
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Appendix 

A. Proof of Lemma^ 

Note that u 2 D2(u) is polynomial with order 2 at least and order 4 at most. For 
k = 1 , 3, let Ck be the coefficients of u in the polynomial u 2 D2(u). Then, 

c\ = -2(1 -poo)/>(-i)o-4(p(-i)(-i)Poi +/>(-i)iA)(-i)) <0, 
C3 = -2(l-poo)pio-4(pi(_i)/7oi+/>iiA)(-i)) <0. 

Hence, if both of u[ l ' m ^ and — Wj max ^ are the solutions of u 2 D2(u) — 0, then 

_/ (l,max) / (l,max)x3N , (l.maxjo/r-, / (l,max)x „ / (l.max)^ „ 

2(au\ '+C3(u\' 'y) = (u\' ') [D 2 {u\ ')-D 2 {-u\ ))=0. 



40 



M. Kobayashi and M. Miyazawa 



Since u\' > 0, this holds true if and only if c\ = C3 = 0, which is equivalent 
to that pq\ = Po(-l) — P(-i)o = P10 = because poo = 1 is impossible. Hence, 
u 2 D2(u) = has the two solutions «j 1,max ' and — Mj 1,max ' if and only if (v-a) does 
not hold. In this case, we have c\ = C3 = 0, which implies that ip-D^iu) is an even 

function. Since u 2 D2(u) = has only real solutions including u^ ,mBX ^ by Lemmas^ 
and [5] we complete the proof. □ 

B. Proof of Lemma^7\ 

By \Y1) , we have 

I I PijH I- I E (-l) i+ ^V = l. 

i£{-i,o,i} ;e{-i,o,i} ie{-i,o,i}M-i,o,i} 
Subtracting both sides of these equations, we have 

pmx + poiy + Poi-i^ 1 +P(-i)o*~ 1 =0. 
Since x,y are positive, this equation holds true if and only if 

P 10 =Poi = Po(-\) =P(-i)o = 0- 
This is the condition that (v-a) does not hold. □ 



